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We discuss the statistics of tunnelling rates in the presence of chaotic classical dynamics. This 
applies to resonance widths in chaotic metastable wells and to tunnelling splittings in chaotic sym- 
metric double wells. The theory is based on using the properties of a semiclassical tunnelling operator 
together with random matrix theory arguments about wave function overlaps. The resulting distri- 
bution depends on the stability of a specific tunnelling orbit and is therefore not universal. However 
it does reduce to the universal Porter-Thomas form as the orbit becomes very unstable. For some 
choices of system parameters there are systematic deviations which we explain in terms of scarring 
of certain real periodic orbits. The theory is tested in a model symmetric double well problem and 
possible experimental realisations are discussed. 



PACS numbers: 03.65.Sq, 73.40Gk, 05.45. Mt, 05.45. 

Tunnelling is crucial in describing many physical phe- 
nomena, from chemical and nuclear reactions to con- 
ductances in mesoscopic devices and ionisation rates in 
atomic systems. When such systems are complex, it is 
natural to model tunnelling effects using random matrix 
theory We show here that when the underlying sys- 
tem is one of clean chaotic dynamics, successful statisti- 
cal modelling demands explicit incorporation of nonuni- 
versal, but simple, dynamical information. We derive a 
distribution for the tunnelling rate which depends on a 
single parameter, calculated from the stability properties 
of the dominant tunnelling orbit. 

The signatures of chaos in tunnelling rates have been 
receiving a growing amount of attention, two important 
regimes having been considered. The first is that the 
quantum state is initially localised in a region where the 
dynamics is largely nonchaotic and one wants to under- 
stand the tunnelling rate through chaotic regions of phase 
space [|]-(|. The second, and the one we shall focus on, 
is that virtually all of the energetically accessible phase 
space is chaotic |5],^| so that the quantum state is initially 
localised in a chaotic region of phase space. These two 
situations are different in many important ways. In par- 
ticular, the statistical distribution of the tunnelling rates 
in the first regime has power law decays || whereas we 
show that the distribution in the second regime has expo- 
nential decay. The result is a generalisation of the Porter- 
Thomas distribution |q] used to model neutron and pro- 
ton resonances |7H9(1 and conductance peak heights in 
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It is shown in Refs. Q that the average tunnelling rate 
from an energetically-connected region of phase space is 
determined by a complex orbit we will call the instanton. 
Fluctuations about this average appear to be pseudo- 
random in the chaotic case and are given by properties 
of the wavefunction in an area localised around a real 
extension of the instanton 101 . To characterise these fluc- 



tuations we define a rescaled tunnelling rate as follows. 
In the case of metastable wells the absolute tunnelling 
rate of a given state labelled by n is measured by the res- 
onance width, or inverse lifetime r„. The corresponding 
normalised tunnelling rate is defined to be 



Vn 



r„/f, 



(i) 



where T(E, ft) = (T n ) is a local average computed for a 
given set of physical parameters. T(E, ft) is a smooth, 
monotonic function of its arguments and is given by an 
explicit formula in terms of the (purely imaginary) ac- 
tion and stability of the instanton j5j . A similar defini- 
tion holds for splittings in double wells and in either case 
(yn) = 1 by construction. 

Fluctuations in y n are calculated using a tunnelling 
operator, 7~, which is constructed from the semiclassical 
Green's function and can be interpreted as transporting 
the wavefunction across the barrier. Specifically, 



y n cx (n\T\n), 



(2) 



where \n) is the wavefunction (which may be calculated 
while ignoring tunnelling effects) represented in a Hilbcrt 
space which quantises a surface of section. A closed-form 
expression can be found by expanding classical actions 
in quadratic order around the instanton. Details can be 
found in Q, but for present purposes it is enough to 
know its spectrum. For a two-dimensional system this 
is {A* |A|^ 2 , k = 0, 1 • • •}, where A is the inverse of the 
stability of the instanton orbit, is always less than unity 
in magnitude and can easily be found using real dynam- 
ics in the inverted potential. (The discussion is readily 
generalised to higher dimension but we refrain from do- 
ing so for clarity.) To derive statistical distributions for 
y n , we will make statistical assumptions about the state 
|n), but not about the tunnelling operator. The resulting 
distribution depends parametrically on A and is therefore 
system-specific and not universal. 
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We can always express T as a diagonal operator in its 
own eigenbasis: J2k ^ k I We then have 



fc=0 



a > A fe |(fc|n}| 2 



k=Q 



k, x k \ 2 



(3) 



where we denote Xk = (k\n) and the prefactor a = 1 — A 
ensures that (y) = 1. The states \n) are normalised so 
that on average \xk\ 2 is unity. We now make the statis- 
tical ansatz that the overlaps can be treated as Gaussian 
random variables. This is the basis of almost all sta- 
tistical treatments of wave functions, going back to the 
seminal work of Porter and Thomas . In the event that 
there is a time reversal symmetry, Xk can be expressed 
as a single real number, leading to GOE statistics. If 
there is no such symmetry then Xk is complex and will 
be described by two statistically independent quantities 
leading to GUE statistics. 

We simplify the derivation by assuming GOE statistics; 
the generalisation to GUE is simple and we give the final 
result for both. We start by assuming that the Xk are sta- 
tistically independent and given by the joint distribution 
P(x)dx = Y[ k [exp(— x 2 ,/2)/s/2tt\ dx k where x = {x k }- 
Wc then note that 



/ dxP(x)^(y- a ^A fc 4 

J k=0 



(4) 



We use the identity S(z) = J dt exp(itz)/27r in the above 
expression and observe that each Xk involves a simple 
Gaussian integral. The final result (and generalising to 
the GUE case) is 



2tt 



OO 

dt e lty J| ( 

k=0 



Oj \ -0/2 

I + -a\H) , (5) 



where (3 — 1 and 2 for GOE and GUE respectively. The 
product above converges rapidly provided A is not too 
close to unity so the formula can easily be used to calcu- 
late P(y; A) in practice. 

This result has a simple interpretation if we think of 
each eigenstate of T as providing a distinct and statisti- 
cally independent channel to tunnel; it then corresponds 
to the discussion of Ref. |12| but with an infinite num- 
ber of distinctly weighted open channels. We could even 
imagine adding a weak magnetic field so as to interpolate 
between the GOE and GUE limits as in jTJj although we 
refrain from that here. As mentioned, we have (y) = 1 
by construction; the second moment is 



(y > = 1 + « 



2 1 - A 



/31 + A 



(6) 



The channel interpretation helps in a qualitative under- 
standing of this distribution as we vary A. For small 
A, only the first (k = 0) channel plays any significant 



role and the distribution is of the Porter-Thomas form: 
exp(— y/2)/y/2ny and exp(— y) for fi = 1 and 2 respec- 
tively. This can be understood analytically from (||) by 
doing a branch-point / residue analysis around the nearest 
singularity at t — i(3/2a. This distribution is commonly 
used to model point tunnelling contacts ]l0| , ^2| , ^3[ . It 
is often a very accurate approximation but its validity 
is not universally guaranteed, as we shall discuss. For 
A close to unity, many channels contribute significantly, 
the fluctuations around the mean are accordingly reduced 
and the distribution approaches a Gaussian with variance 
<7 2 ~ a I '(3 (which becomes a delta function for small a). 

For A > and y < we close the contour of (||) in the 
lower half plane; since the integrand has no singularities 
there, the result is simply zero. This is consistent with 
the fact that T is a positive definite operator so that 
Eq.(||) does not admit the possibility of negative y. By 
the same argument, any derivative of p(y) is also zero 
for y < implying a nonanalyticity at y = with p(y) 
going to zero faster than any power of y for y small and 
positive. In the opposite limit j » ) we can expand 
around the first singularity to obtain 



P(y, A)goe 
P{y, A)gue 



exp(— y/2a) 

y/2nay 
exp(-y/a) 



(7) 



This falls off exponentially with y, and not with a power 
law as observed in the chaos-assisted regime [|| . 

Equations (^|) , (|^) and (^|) remain valid when A is neg- 
ative. This situation arises when we compute splittings 
in symmetric double wells for which the symmetry is in- 
version through a point rather than reflection through an 
axis Q . In this situation T is no longer positive definite 
and we admit the possibility of negative splittings (for 
which the odd member of a doublet has a lower energy 
than the even member). The distribution then allows all 
values of y, positive and negative. It decays exponentially 
for |y| ^ A as in Eq. (Q) but with different exponents for 
y > and y < because on doing the integral (||) we 
must switch from closing the integration contour in the 
lower half plane to the upper half plane as y changes sign. 
For a given state, we can typically induce a zero splitting 
by tuning one system parameter. A zero splitting means 
that we can construct states which remain localised in 
either well for all time as in the one dimensional time- 
dependent system considered in jTJ|] . 

At this point we contrast our results with standard ran- 
dom matrix theory modelling. In section VII. H of their 
extensive review |8j , Brody et al. show under rather gen- 
eral assumptions that one expects a Gaussian distribu- 
tion for the expectation values of an arbitrary operator 
by showing that all of the moments of the distribution ap- 
proach those of a Gaussian. This can be understood as a 
sort of central limit theorem. One of their assumptions is 
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that the operator is non-singular, i.e. does not have many 
zero eigenvalues. Because of the exponential decay of the 
eigenvalues of the tunnelling operator T, it is effectively 
singular. Therefore their conclusions do not apply to our 
situation and we have non-Gaussian distributions, ft is 
interesting to note that in the limit |A| — > 1, the operator 
T has an ever increasing number of significant eigenval- 
ues and the distributions do in fact approach Gaussians, 
in conformity with their general considerations. 



of parameters such that the classical dynamics is almost 
fully chaotic. We also show the distribution (||) with 
(3 = 1 and using the corresponding values of A. Clearly, 
the numerically computed histograms are well captured 
by the theoretical distribution. Wc show, for compari- 
son, the Porter-Thomas distribution which clearly fails 
to correctly model the numerical data. We remark that 
this sort of agreement was observed for most parameter 
values as long as the dynamics was fully chaotic. 
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FIG. 1. Results for two typical potentials, with 
(fj,,v,(r) = (0.15,0.17,0.00) above and (0.25,0.50,0.00) be- 
low. In both cases the energy is lower than saddle maximum 
by an amount 0.1. The continuous curves are the theoret- 
ical distributions calculated using the appropriate values of 
A (shown). The dashed curves show the Porter-Thomas dis- 
tribution for comparison. The insets show the corresponding 
instanton orbits and their real extensions. 

Since it is a simpler numerical task to calculate many 
splittings in a double well than to calculate many reso- 
nance widths in a metastable well, we use the former to 
test our predictions and note that any conclusions apply 
identically to the latter. Consider the potential 

V(x, y) = {x 2 - l) 4 + x 2 y 2 + m) 2 + vy + ax 2 y. (8) 

There is a reflection symmetry in x and if the energy is 
less than 1 — i/ 2 /4/i the motion is classically confined ei- 
ther to x < or to x > 0. It is convenient to work at 
fixed energy in order to keep A constant and we do this by 
quantising q — 1/h ||, that is, by finding those values of 
h which are consistent with a specified choice of param- 
eters and energy. This is effectively what happens, for 
example, in scaling problems such as a hydrogen atom in 
a magnetic and electric field |l5| . In Fig. |l| we show his- 
tograms constructed from the q-spectra for two choices 
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FIG. 2. As in Fig. 1 but for parameter values 
(/i, v, a) = (0.25, 0.40, 0.254) for which the real extension of 
the instanton orbit is periodic and the resulting distribution 
is significantly different from the RMT prediction. 

In Fig. |^ we show an exception to the general agree- 
ment. In this case the numerical histogram is interme- 
diate between the theoretical distribution (||) and the 
Porter-Thomas form. We attribute this to the effects 
of scarring [jl6 17|, as follows. The instanton has real 
turning points where the momentum vanishes and the 
position is real. At these points we can elect either to 
integrate in imaginary time in which case the instanton 
retraces itself or to integrate in real time, in which case 
we get a real trajectory. We refer to this real trajec- 
tory as the real extension of the instanton. There is 
no reason why this real extension should itself be pe- 
riodic. Typically it is not. However, the parameters 
of Fig. H have been tuned so that the real extension is 
in fact periodic. We find in this case that the overlaps 
Xk = (k\n) are no longer distributed according to the 
Gaussian P(xk) = exp(— x 2 ./2)/\/2tt as assumed in our 
derivation — there are relatively more large overlaps and 
more small overlaps. This effect can be explained using 
a recent theory of scarring |17| which describes how the 
overlaps between a wavepacket placed on a periodic orbit 
and the chaotic eigenstates deviate from random matrix 
theory. In our problem the eigenvectors \k) behave like 
wavepackets of this type when the real extension is pe- 
riodic. The effect of this deviation from random matrix 
theory is to give more large splittings and more small 
splittings than (|^) predicts and, therefore, to push the 
distribution in the direction of the Porter-Thomas form. 
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We remark that for v = a = 0, the real extension is al- 
ways a periodic orbit and we see anomalous statistics 
for this situation as well. 
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discernible even for rather small values of A. 

Another possibility for nonuniversal statistics would 
be a situation analogous to Fig. 0, where the tunnelling 
route is directly connected to a real periodic orbit. This 
geometry could be engineered into quantum dots and is 
present in the hydrogen atom problem |jl5). In this case 
we predict deviation from the predictions of random ma- 
trix theory. This would be similar in spirit to the recent 
work of Narimanov et al. jl9| who look for dynamical ef- 
fects in the correlations of conductance peaks of quantum 
dots. 



FIG. 3. Results for the inversion-symmetric potential with 
(fi, t) = (0.1, 1.0) for which zero and negative splittings are 
allowed. 

The final case we discuss is if the term vy + ax 2 y in (||) 
is replaced by rxy. Now the potential is symmetric un- 
der (x, y) — ► {—x, —y) rather than under (x, y) — > (—x, y) 
(this symmetry would persist if we were to add a uniform 
magnetic field). In this case A < || and negative split- 
tings can occur. We present the results for a typical case 
in Fig. |[ Again, the theoretical distribution agrees with 
the histogram. 

Our results are relevant to situations in which particles 
tunnel out of or between chaotic regions separated by an 
energetic barrier. Applications include hydrogen atoms 
in parallel electric and magnetic fields where the compe- 
tition between the imposed fields and the Coulomb force 
causes chaotic motion while the presence of the electric 
field causes tunnelling |u| . Dissociative decays of excited 
nuclei and molecules may also fall into this regime. 

Another application is to conductances of quantum 
dots. In the Coulomb blockade regime electrons must 
tunnel into and out of dots which are thought to be 
chaotic. Such experiments have been done |n| leading 
to results which are consistent with the Porter-Thomas 
distribution for the tunnelling widths, just as for the neu- 
tron and proton resonance widths. We contend that the 
reason is that the instanton path in all cases is very un- 
stable, leading to a small value of A. For energies near 
the saddle A ~ exp(— 2 r KLO y /uj x ) [jbSf where ui x and uj y are 
the curvatures of the potential saddle along and trans- 
verse to the instanton, respectively. This is often small, 
but by making the saddle flat in the transverse direction 
or sharp in the instanton direction, it is possible to have 
A be of order unity. It is an interesting question whether 
this can be arranged for the quantum dots. One feature 
which helps in this regard is that we predict a distribution 
which vanishes for small y whereas the Porter-Thomas 
distribution diverges as l/^/y- This difference could be 
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